﻿ Sampling •time-domain discretization •Sampling of bandlimited signals •Reconstruction of analog signals from the sampled signal 1 Sampling theorem 1⎡Δ Δ⎤ ⎛⎞⎛⎞ ut t t=σ+−σ− () ⎜⎟⎜⎟ Δ ⎢⎥ 22Δ ⎝⎠⎝⎠ ⎣⎦ A sample of x(t) is obtained multiplying the continuous-time signal x(t) with the rectangular impulse u(t ) of area 1: Δ 0xtu t x u t≅ ()()()() ΔΔ Another sample can be obtained using the same rectangular impulse shifted by kT s −≅ − ()()() () sssxt u t kT x kT u t kTΔΔ 2 1 sampling process : ∞∞ −≅ − ()() ()() ∑∑ sssxtutkT xkTutkTΔΔ ∞ =−∞ kk=− For small values of Δ: Dirac distribution limut t=δ ()() Δ→ 0Δ ∞∞ limu t kT t kT x kT t−=δ− =δ () ()() () ∑∑ 0sTsssΔ Δ→ kk=− ∞ =−∞ Periodic Dirac distribution ideal sampling of the signal x(t) : ∞ ˆxtxt t xkT tkT=δ = δ− () () ()()( ) sTss∑ ∞ k=− 3 Ideal Sampling Mathematical model ∞ xtxt t xkT tkT=δ = δ− () () ()()( ) sTss∑ ∞ k=− System model 4 2 Spectrum of Ideal Sampled Signal ideal sampled signal: ∞ xtxt t xkT tkT=δ = δ− () () ()()( ) sTss∑ ∞ k=− the spectrum of the ideal sampled signal : ˆ Xxt tω= δF ( ) () () {} sT 5 The Fourier transform of a product is convolution of the Fourier transforms The Fourier transform of the periodic Dirac’s distribution is also a periodic Dirac’s distribution ∞ 12 2⎛⎞ππ ˆ δ = ω∗ δω− XxttX kω=F ( ) () ()() {} ∑ sT⎜⎟ 2TTπ =−∞ kss⎝⎠ ∞∞ ⎛⎞ 1212⎛⎞ππ ˆ XXk Xkω= ω∗δω− = ω− ()() ⎜⎟ ⎜⎟ ∑∑ ∞ =−∞ kksss sTTTT=− ⎝⎠ ⎝⎠ The effect of the convolution of a specified function with the periodic Dirac’s distribution is the periodic repetition of the considered function 6 3 ∞ − ()()( ) ∑ ssxtxkTtkT=δ ∞ k=− ∞ 12⎛⎞π ˆ XXkω= ω− () ⎜⎟ ∑ ∞ kssTT=− ⎝⎠ The spectrum of an ideal sampled signal is the periodic repetition of the spectrum of the original signal The period is inverse proportional with the sampling step T s 7 Spectrum of original signal Spectrum of periodic Dirac distribution Spectrum of ideal sampled signal Aliasing ∞ 12⎛⎞ π ˆ XXkω= ω− 8error() ⎜⎟ ∑ ∞ kssTT=− ⎝⎠ 4 Sampling Band-limited Signals x(t)-band-limited () 0 if ≥ω=ω MXω 9 • successive replica of the original signal spectrum are not superimposed and the spectrum of the original signal can be recovered by low-pass ideal filtering for 2ω≥ω sM 10 5 Ideal Low-pass Filter tsinω c () () () thpH= ↔ω=ω cω tπ 11 the aliasing error can be avoided Sampling freq 2ω≥ ω sM Cutoff freq for the low pass filter ω≤ω−ω csMMω≤ For perfect reconstruction : 0HT= () rs 12 6 The frequency response of the reconstruction filter is: T⎧ω≤ω, cs⎪ ≤ω ≤ω −ω ω= ω= ()() ⎨ csMMω rsHTpω cω0, ⎪ c>ω ⎩ Its response : ˆxtxtht=∗ ()()() rr with the spectrum: ˆ ω⋅ ω ()()() rrXXHω= ∞ 1 XkTp X= ω− ω ω = ω ()() () ∑ cssω ksT=− ∞ , xtxtaew= ()() r 13 ω When the condition: is not verified, the aliasing error2Msω≥ appears ω−ω ω sxnT n∈ M), it is uniquely determined by its samples { if the sampling frequency is higher or equal than twice the maximum frequency of the signal: 2ω≥ω sM the original signal can be reconstructed from its samples a e w: sin2tkT∞ω− () ccsω ()() ∑ rsxt xkT= ωω− () ∞ kscstkT=− if the cut-off frequency of ideal low-pass reconstruction filter : ω≤ω ≤ω −ω 17csMM Ideal Low-pass Filtering Reconstruction The signal reconstructed from curves of type sin x / x Interpolation 18 9 • In the sampling moments, only one sample defines the signal sinTnk∞ω− () Ms () () ∑ ssxnT xkT= ω− () kMsTnk=− ∞ sinnk∞∞π− () sω ⇒ = = δ = xnT x kT x kT x nTω= () ()() () ,ssMssnk∑∑ 2kknk=− π− () ∞ =−∞ 1, for n k= ⎧ δ= ,nk⎨ 0, for n k≠ ⎩ • In between the sampling points, infinite number of terms 19 Linear interpolation • approximates the signal using straight lines that unify points determined by the samples • Reconstruction filter is triangular • errors 20 10 2 2 ω ⎛⎞ ⎛⎞ sinsTω sinπ ⎜⎟ ⎜⎟ 2sω ω= = () ⎜⎟ rs sHT T⎜⎟ ⎜⎟ sTωω ⎜⎟ π ⎝⎠ 2⎜⎟ sω ⎝⎠ ideal 21 Reconstruction by Zero Order Extrapolation TTωω 2sin sinTTssssω jjT−−ω 2222s⎛⎞ ⎯→= () Trsht p t e e T=− s← T⎜⎟ ω ⎝⎠ 22sω 2 ω ω sinsTω sinTsπ j−π j−ω ω 22ssω HeT eω= () rs = sTωω 22π 2ω s 11 • Difference between filters is significant ⇒ large errors The frequency response of the reconstruction filter is not flat in the pass band 23perfectly The spectrum of the reconstructed signal: ω sinωπ j∞ −π ω ssω ω ω= ω−ω () () () () ∑ rrsXXH e Xkω= ω k=− ∞ π sω 24 12 The reconstructed signal is not band-limited The reconstruction error can be further reduced if an ideal low-pass filter is added ω Even after low pass filtering errors persist It is decreased if ω> sM> 25 Ideal Sampling of Periodic Signals • periodic signal of periodT, band limited at theN’th harmonic (maximum frequency) 2π Nω=; ω ω= 00M T 0 • Ideal sampling with sampling frequency MMω= ω ∈, 0s • Effect: periodic repetition of its spectrum, weighted by inverse of the sampling step • Spectral components should not overlap: NNMNω ω ⎩ NNω ω 00s 22Nω>ω=ω 0sM Reminder: sampling non-periodic signals 2ω sM≥ω On the period of the fastest harmonic component more than 2 samples (at least 3) must be taken 29 T- period of the fundamental 0 ideal sampling with ω s=(2N+R)ω0then: 22ππ 2NR=+ () 0sTT 0T ; 1, 2,TR== s… 2NR+ 2N+R = maximum different samples on a single fundamental period of the original signal 30 15 Principle of Sampling Oscilloscope Same result : successive samples from successive periods = + () ()() 00sssxkT x T kT x KT kT=+ 1 ′ =+ 000ssTKTkTKT T=+ 2NR+ 31 2Nω>ω http://www jhu edu/~signals/sampling0s NNω > Mω tΔ ⇒ Minimum duration of impulses, Δt reconstruction by zero order extrapolation : particular case of PAM, for Δt=Ts 36 18 Natural sampling Used in time division multiplexing (TDM) of analog signals switch driven by a periodic rectangular train of impulses high impulse : input signal is transmitted 37 low impulse : zero xt xtq t xt ht t⎡⎤ ==∗δ () () () () () () ssTT⎣⎦ ∞∞ xthtkT xthtkT=− () ( ) () ( ) ∑∑ ss=− ∞ =−∞ kk=− •input signal multiplied with the periodic impulse train q(t) Ts tωΔ 2sintωΔ jt− ⎛⎞ 22Δ =−←⎯→ω= ()() ⎜⎟ tht p t H eΔ ⎝⎠ 22ω 38 19 Natural Sampled Signal Spectrum Xxthtt⎡⎤F ω= ∗δ ( ) () () () {} sT⎣⎦ tωΔ ⎡⎤ sintωΔ − j∞ ⎢⎥ 2122π Xet k=ω Δ ⋅ δω−ω ()() ∑ ⎢⎥ s∗ ωΔ 2ktT=− π s∞ ⎢⎥ 2⎣ ⎦ ktωΔ sinsktωΔ s∞ jt− 22Δ Xe k=ω ⋅δω−ω () () ∑ s∗ ωΔ skktT=− s∞ 2 numerical factor ktωΔ sinsktωΔ s∞ jt− 22Δ e= −ω () ∑ sXk⋅ω sktTωΔ ∞ sk=− 39 2 numerical factor sT⎧ ≤ω ⎪ c,ω Htω= Δ () ⎨ •No distortion of central spectral lobe, as 0,⎪ω>ω c⎩ opposed to sampling and hold •Reconstruction by ideal low-pass filtering like ≤ω≤ω−ω in the case of ideal sampling McsMω 40 20 Connection between spectra of discrete and corresponding continuous signals xt X xn xnT X←⎯→ω = ←⎯→Ω; ()()()() [] aa asdd ∞ 12π ˆ xt xt t X X k=δˆ, ←⎯→ω= ω−ω ω= () () () () ( ) sassT∑ ∞ kssTT=− ∞ ⎫ ˆ⎧ xt xt t X x t t kT=δˆ ←⎯→ω= δ− F () () () () ( ) () ⎬ ∑ sasT⎨ k=− ∞ ⎩⎭ ∞∞ ⎧⎫ xt tkT xkT tkT=δ F = δ− ()() ()() ⎨⎬ ∑∑ as ass− kk=− ∞ =−∞ ⎩⎭ ∞ xkT t kT=δ F () ( ) {} ∑ as s− k=− ∞ ∞ 41sjkT−ω () xkTe= as∑ k=− ∞ ∞∞ 1j ˆskT−ω ω− ω = () () () ∑∑ as asXXkxkTeω= ∞ =−∞ kksT=− The Discrete Fourier spectrum: ∞∞ jk jk−Ω−Ω = () () [] ∑∑ asddXxkexkTeΩ= kk=−∞ =−∞ It can be observed that: ∞∞ − Ω sjkT jk−ω Ω ()() ∑∑ as asx kTe x kTe= ω= ∞ =−∞ kk=− sT Connection between spectra of discrete and corresponding continuous signals : ∞ 12⎛⎞Ωπ − () ⎜⎟ ∑ adXXkΩ= ∞ ⎝⎠ 42ksssTTT=− 21 The link between the two frequency axes of spectra of the discrete signal and of the corresponding continuous signal is: =ω sTΩ ⇒spectrum periodic() dXΩ The maximum frequencies are also related: π TTΩ=ω ≤; 43ssMMω M Sampling of discrete-time signals Unknown maximum frequency of a band-limited signal generally ⇒ use maximum sampling frequency allowed Problem: the sampling frequency used is too high (oversampling) Solution: reduce the number of samples after the signal was sampled = Sampling a discrete-time signal ∞ ˆxnxn nxn nkN=δ δ− [] [] [] [] [ ] ∑ N = ∞ k=− ∞ xkN n kN=δ − [][ ] ∑ ∞ k=− 44 22 sampling step N=3 45 Spectrum of Discrete-time Sampled Signal ∞ 2π ⎯→Ω δ Ω = Ω δ Ω− Ω Ω = nkδ←; ()() [] ∑ sNss ssΩ kN=−∞ ˆ δF () [][] {} NXxnnΩ= ∞ 12 2π π Xk=Ω ; δΩ−ΩΩ= ()() ∑ ss⊗ 2NNπ =−∞ k Spectrum of a discrete-time sampled signal = periodic repetition of the spectrum of the original signal 11NN−− 112π ˆ XX k XkΩ= Ω∗ δΩ−Ω = Ω−Ω Ω=, ()()()() rr rsss∑∑ 00kkNNN== 46 Restriction for one period, 2π 23 sampling step N=3 Spectrum periodic of period Ω s 47 The continuous-time signal, x(t), band-limited at ω, ideal sampled, M sampling step T⇒ sband-limited discrete-time signal x[n], also sampled, step N The spectral lobes do not overlap if the analog signal x(t) could have been sampled with the sampling step T’=NTfollowing the sampling ss theorem Initially the signal was oversampled /N, the perfect Because the maximum frequency isΩ M=ωMTs, ΩS=2π reconstruction condition is: 2Ω Ω s≥M ππ ≤≤= NTTTNT′′; 48ssss MMωω 24 Aliasing error Ω−Ω M-1) ⇒ periodic signal , period[ sπΩ Reconstruction condition:If , no overlapping MN ≥ •multiplication with the temporal window: 2π ⎧ , 0 1nN≤ ≤− ⎪ x n xn xnw n w nN==; [] [] [] [] [] ⎨ rrr = ⎪0, otherwise ⎩ Low sampling step, N<M ⇒ successive temporal groups overlap; aliasing error appears Original signal can not be reconstructed using a temporal window 58 29 Practical Sampling of Continuous-time Signals •bandwidth of the signal to be sampled not known •high frequency spectral components, not important for the application (e g produced by the noise that perturbs the useful signal) •aliasing errors Avoid using low-pass filter: anti-aliasing filter before the sampling system 59 •digital telephony system (audio signal) : •maximum frequency f=3 4 kHz M •Theoretical sampling frequency f=6 8 kHz s •sampling frequency f=8 kHz s •television system (video signal): •maximum frequency f=5 MHz M •Theoretical sampling frequency f=10 MHz s •Sampling frequency f=18 MHz s 60 30 originalspeech wav Speech signal without aliasing aliasingspeech wav Speech signal with aliasing music wav Music signal without aliasing Music signal with aliasing aliasingmusic wav 61 Sampling of Band-pass Signals Spectrum of low-pass signals centered on 0 Spectrum of band-pass signals has a support of the form supp X ,ω=−ω −ω ∪ω ω , ⎤⎡ ⎤ () {} mmMM⎡ ⎣⎦⎣ ⎦ Perfect reconstruction of ideal sampled band-pass signal using WKS theorem: 2ω ω s≥M Occasionally, perfect reconstruction of ideal sampled band-pass signals can be achieved even if a sampling frequency inferior to the Nyquist frequencywas 62used 31 Narrow-band Band-pass Signals −ω mMω Narrow-band band-pass signals : 1< mω Support of spectrum for ideal sampled narrow-band band-pass signal supp X , ,nn nnω =∪ −ω +ω−ω+ω ∪ω+ωω +ω ⎡⎤ () {} {} mmMMsss ss⎡⎤ ⎣⎦⎣⎦ n∈ Here: ω= ω sM Perfect reconstruction: ideal band-pass filter even for sampling 63 frequency inferior to the Nyquist frequency Ideal Band-pass Filter ⎧ω<ω<ω1, 12cc⎪ = ω− ω Hppω= ()() () ⎨ 21ccBPωω 0, otherwise⎪ ⎩ ω 21sin sinccttω ht=− () BP ttππ 64 32 Perfect reconstruction condition : ,, ,kk ll kl−ω + ω −ω + ω ∩ ω + ω ω + ω = ∅ ∀ ∈ ⎡⎤ mmMMss ss⎡⎤ ⎣⎦⎣⎦ Translation of interval Translation of interval [-ω m,ωM] with lωs M,-ωm] with kωs[ω For l=0: kk k−ω+ω−ω+ω∩ωω =∅∀∈, , , ⎡⎤ mmMMss⎡⎤ ⎣⎦⎣⎦ ω ≤ ω ⎧ 22mk−ω ω ω MmMs + ω≤ ⎨ s⇔≤ + ω ≥ ω 11kkk−ω + () MMs + ⎩ If there are integer values of k for which this condition is satisfied then there are values for the sampling frequency inferior to the Nyquist frequency, for which the narrow-band band-pass signals can be reconstructed from their samples 65 mω ≤= 00kn< ω mMω− sampling frequency must belong to intervals of the form: 22ω ⎡⎤ mMω … {} with 1, ,kn∈0, ⎢⎥ 1kk+ ⎣⎦ Example ω 4mn== 10ω=π π If and then ⇒ k =1,2,3, 4 8mω= M0ω− mMω The corresponding intervals for the sampling frequency are: ω∈ π∪ π π∪ π π∪ π π∪ π ∞4 5 , 5 33 6 66 , 8 10 , 16 20 , {} [][][][] s 66 33